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AND FAGNER BERNARDINI RODRIGUES 


Abstract. We consider stochastic processes arising from dynamical systems by evaluating 
an observable function along the orbits of the system. The novelty is that we will consider 
observables achieving a global maximum value (possible infinite) at multiple points with spe¬ 
cial emphasis for the case where these maximal points are correlated or bound by belonging 
to the same orbit of a certain chosen point. These multiple correlated maxima can be seen 
as a new mechanism creating clustering. We recall that clustering was intimately connected 
with periodicity when the maximum was achieved at a single point. We will study this 
mechanism for creating clustering and will address the existence of limiting Extreme Value 
Laws, the repercussions on the value of the Extremal Index, the impact on the limit of Rare 
Events Points Processes, the influence on clustering patterns and the competition of domains 
of attraction. We also consider briefly and for comparison purposes multiple uncorrelated 
maxima. The systems considered include expanding maps of the interval such as Rychlik 
maps but also maps with an indifferent fixed point such as Manneville-Pommeau maps. 

1. Introduction 

In the past few years, the study of Extreme Value Theory (EVT) for dynamical systems 
has been a subject of much interest. We mention in particular the inspiring paper of Collet 
fColoT] and refer to |Frel3| for a review on further developments and references. Extreme 
events (which occur with small probability and for that reason are also called rare events) 
are characterised by abnormally high observations that are identified as exceedances of high 
thresholds. 

In this dynamical setting the observation data comes from starting the system at a certain 
initial state and evaluate a certain observable function along the subsequent states through 
which the system goes while time goes by. Then, exceedances of a high threshold correspond 
to entrances or hits of the orbits of the system to some designated target sets on the phase 
space. As the threshold increases, the respective target sets shrink. This fact explains the 
connection between the existence of distributional limits, called Extreme Value Laws (EVL), 
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for the partial maximum of observations when the thresholds increase to the possible maximum 
value, and the existence of distributional limits, called Hitting Times Statistics (HTS), for the 
waiting time before hitting the corresponding target sets as they shrink. This connection 
hinted in |Col01| was formally established in |FFT101 IFFT11| . 

In the study of HTS, in most cases, the neighbourhoods are either cylinder sets or metric 
balls and shrink to a point £ in the phase space. In the study of EVL the observable 99 has 
a global maximum at £ and, in almost all cases, the exceedances correspond to metric balls 
around £. The limiting laws for HTS and EVL were proved to be the same in |FFT101 IFFTllj 
and typically one obtains a standard exponential distribution H{t) = 1 — e -r , with r > 0. 
Typically, here, is used in the sense that in most results it is shown that for almost every 
point £ (with respect to the invariant measure) one gets a standard exponential HTS and 
EVL. See for example |HSV99I ICdOll IBSTV0.11 IHNT121 ICClll lPSl4] and |Sau09] for an 
excellent review. 

Following the work of Hirata |Hir91| on Axiom A diffeomorphisms, it is known that at 
periodic points, i.e., when £ is a periodic point, a parameter 0 < 9 < 1 appears in the 
asymptotic distribution: H(t) = 1 —e _0r . This observation was further developed in the paper 
pwno] . where a compound Poisson distribution for the number of hits to target sets composed 
of unions of dynamical cylinders was obtained for i[) mixing measures. Then in |FFT12j . using 
the relation between HTS and EVL, a new technique based on suitably adjusted dependent 
conditions allowed to show the existence of EVL (and consequently HTS) in the case £ is a 
periodic point and the target sets are metric balls for systems with a strong form of decay 
of correlations. Due to the connection with EVL the parameter 0 < 9 < 1 appearing in the 
limit distribution H(r ) = 1 — e -0T was identified to be the Extremal Index (El), a parameter 
already appearing in the context of classical EVT and which can be seen as the inverse of the 
average cluster size. 

Moreover, in |FFT12j . for uniformly expanding systems such as the doubling map, a di¬ 
chotomy was shown which states that either £ is periodic and we have an El with a very precise 
formula depending on the expansion rate at £, or for every non-periodic £, we have an El equal 
to 1 (which means no clustering). The dichotomy was then obtained for more general sys¬ 
tems such as: conformal repellers |FP12] . systems with spectral gaps for the Perron-Frobenius 
operator |Kell2| and systems with strong decay of correlations |IAFV14j . 

The study of EVL and HTS can be enhanced by considering Rare Events Point Processes 
(REPP). These point processes keep record of the exceedances of the high thresholds by 
counting the number of such exceedances on a rescaled time interval. At typical points it was 
known that the REPP converge in distribution to a standard Poisson process. From |HV09j . at 
periodic points, we expect that the REPP converge to compound Poisson process instead. In 
|FFT13| . a study regarding the convergence of REPP was performed for metric balls as target 
sets around periodic points £ and for non-uniformly hyperbolic systems. The limit process 
obtained for the convergence in distribution of the REPP in |FFT13j was a compound Poisson 
process which could be described as combination of Poisson process ruling the positions of the 
clusters in the time line being that these positions are marked by the cluster sizes ruled by a 
geometric multiplicity distribution. Both components are determined by the value of El 9. 

Again, as proved in (AFV14) . for continuous systems with strong decay of correlations 
a dichotomy regarding the convergence of REPP holds: either it converges to a compound 
Poisson distribution at repelling periodic points £, with a geometric multiplicity distribution, 
or it converges to a standard Poisson process at every other non-periodic point £. Moreover, as 
already seen in |FFT13j . the clustering pattern observed at repelling periodic points £ obeys 
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a very rigid pattern: it consists of ‘bulk’ of strictly decreasing exceedances observed at precise 
fixed times, corresponding to the period. 

Contrary to the customary case studies, in which the set of maximising points A4 of the 
observable tp is reduced to a single point (, in this paper we will consider the case of multiple 
maxima, with special emphasis for the case of correlated maxima, where J\4 consists of finite 
number of points bound together by belonging to the same orbit. This binding gives a mech¬ 
anism of creating clustering of exceedances, which we will study and explore below. We will 
also consider the case of uncorrelated maxima but the focus will be turned to the correlated 
case. 

We remark that in the literature not all examples regarding the study of rare events are 
reduced to the cases in which the target sets shrink to a point. For example, in |CCC09| the 
authors consider a sets generated by a proper subshift of finite type of a one-sided irreducible 
and aperiodic shift of finite type and in [KL09| the authors consider target sets shrinking 
to the diagonal of the phase space of a system obtained by coupling two expanding interval 
maps. The setting which probably most resembles ours is that studied in |HNT12j . where the 
authors, in the EVL approach, consider observables with multiple maxima that are chosen 
independently as typical points for the invariant probability measure. This way, for these 
uncorrelated maxima, they obtained the standard exponential distribution as limiting EVL, 
which means that there is no clustering of exceedances. 

In classical Extreme Value Theory, it is usual to use normalising thresholds (t! n )n.eN tli a t are 
linear sequences depending on a parameter y, say u n = y/a n +b n . As explained in |FFT10| . the 
behaviour of the observable ip as a function of the distance to ( is responsible for determining 
the type of EVL that applies: Gumbel, Frechet or Weibull. Hence, an interesting question 
that immediately arises when considering multiple maxima is the competition between the 
different types. 

By taking multiple maxima over the same orbit we will see the following consequences: 

(1) appearance of clustering not caused necessarily by periodic orbits; 

(2) the possibility of creating different clustering patterns (not reduced to a bulk of strictly 
decreasing observations over the threshold); 

(3) the possibility of affecting the El when we already consider periodic points; 

(4) different multiplicity distributions for the limit of REPP; 

(5) competition between different types of distributions. 

The systems considered include systems with decay of correlations against L 1 observables 
such as Rychlik maps |Ryc83 [ and piecewise uniformly expanding maps in higher dimensions 
studied in |Sau00j but also non-uniformly expanding maps such as intermittent maps. 

We believe that this study will be a precursor of further developments regarding more 
sophisticated maximal sets A4, that could be, for example, submanifolds of codimension 1. In 
fact, some work in that direction is already being done by some of the authors and we believe 
it carries a large potential of applications since it will more easily accommodate physical 
observables in applications like to meteorology, where the set A4 plays the role of a critical 
set, which is the source of abnormal and possibly catastrophic events. Understanding the 
geometry and the recurrence properties of J\4, as we will see here in a much simpler context, 
will provide knowledge about the extremal behaviour like for example the clustering patterns 
that could help devise early warning systems. 

An example of possible application is to structural failures. The anticipated study of the 
type of the clustering patterns of a certain natural phenomenon is of crucial importance, on 
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one side to predict the likelihood of such failure and on, another side, to help designing the 
material and structures to stand stronger against the natural causes they have to face. 

The unfolding of such possibilities will give deeper understanding of the possible outputs 
of EVT for dynamical systems, as it will open a door for modelling physical phenomena with 
some underlying periodic effect, which is often difficult to perceive. Moreover, again, these 
short recurrence mechanisms have an enormous potential as a source of examples for the 
classical EVT of stochastic processes and serving as a model for several practical situations. 


2. The setting and background 

Take a system (df, B, fi, /), where X is a Riemannian manifold, B is the Borel u-algebra, 
/ : X —> X is a measurable map and p an /-invariant probability measure. Suppose that 
the time series Vo, X \.... arises from such a system simply by evaluating a given observable 
ip : X — > MU {± 00 } along the orbits of the system, or in other words, the time evolution given 
by successive iterations by /: 

X n = p o f n , for each n E N. (1) 

Clearly, Xq, X\, ... defined in this way is not necessarily an independent sequence. However, 
/-invariance of p guarantees that this stochastic process is stationary. 


We suppose that the r.v. p : X 


U {± 00 } has N global maxima £ 1 , • • • ,£/v £ X (we 


allow y>(£i) = + 00 ) and assume they all belong to the orbit of the point /. (See (18) below). 
Suppose <p and p are sufficiently regular in the following sense: 

(Rl) for u sufficiently close to up := p(£i) (i E {1,..., IV}), 


U(u) := {x € X : <p(x) > uj = {Vo > u} 

corresponds to a union of balls centered at the points i.e., U(u) = Ufli B £i (/?;) with 
S{ = Ei(u). Moreover, the quantity p(U(u)), as a function of u, varies continuously on 
a neighborhood of up. 


2.1. Extreme Value Laws. In this paper, we will use an extreme value approach rather 
than an hitting times approach, which we have already mentioned to be two sides of the same 
coin as can be fully appreciated in [FFTlOl IFFT11] . 

We are interested in studying the extremal behaviour of the stochastic process Vo, Vi, ... 
which is tied to the occurrence of exceedances of high levels u. The occurrence of an exceedance 
at time j E No means that the event {Xj > u} occurs, where u is close to uf • Observe that a 
realisation of the stochastic process Vo, Vi,... is achieved if we pick, at random and according 
to the measure p, a point x E X, compute its orbit and evaluate <p along it. Then, saying 
that an exceedance occurs at time j means that the orbit of the point x hits one of the balls 
in U(u) at time j, i.e., f 3 {x) E B Si (£i) for some i E {1,..., V}. 

Given a stochastic process Vo, Vi,... we define a new sequence of random variables M\, M 2 ,... 
given by 

M n = max{V 0 , Vi,..., V n _i}. (2) 

On the independent context, i.e., when the stochastic process Vo, Vi,... is a sequence of 
independent and identically distributed (i.i.d.) r.v., the first statement regarding M n is that 
M n converges almost surely (a.s.) to uf- Then, the next natural question is whether we can 
find a distributional limit for M n . when conveniently normalised. 
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Definition 2.1. We say that we have an Extreme Value Law (EVL) for M n if there is a 
non-degenerate d.f. H : M —> [0,1] with H( 0) =0 and, for every r > 0, there exists a sequence 
of levels u n = u n (r), n = 1,2,..., such that 

n//(Vo > u n ) — y t, as n —> oo, (3) 

and for which the following holds: 

/i(M„ < u n ) -A H(t ), as n —> oo. (4) 

where the convergence is meant at the continuity points of H(t). 

The motivation for using a normalising sequence (un.)neN satisfying ([3]) comes from the case 
when Xq,X\, ... are independent and identically distributed (i.i.d.). In this setting, it is clear 
that n{M n < it) = ( F{u )) n , where F is the d.f. of Xq. Hence, condition ([3]) implies that 

n(M n < u n ) = (1 - n(X 0 > u n )) n ~ p - -J Ae r , (5) 

as n—> oo. Moreover, the reciprocal is also true (see |LLR83l Theorem 1.5.1] for more details). 
Note that in this case H{t) = 1 — e~ T is the standard exponential d.f. 


2.2. The existence of Extreme Value Laws. In |FFT15| . the authors synthesised the 
conditions in |FF08| . in the absence of clustering, and in |FFT12j . in the presence of clustering, 
to a couple of general conditions that apply to to general stationary stochastic processes, both 
in the presence and absence of clustering, which allow to prove the existence of EVL. Moreover, 
these conditions are particularly tailored to the application of dynamical systems and follow 
from a strong form of decay of correlations (against L 1 observables) to be defined below. 

In what follows for every A 6 £>, we denote the complement of A as A c := X \ A. 

For some u 6 M, q € N, we define the events: 


U(u) := {V 0 > u}, 

<? 

Mu) ■= U(u) n p| r\U(u) c ) = {Vo > u, Vi V, < u}. (6) 

2=1 


where A q (u) corresponds to the case where we have an extreme event at time zero that is not 
followed by another one up to time t = q. This is a condition clearly pointing to the absence 
of clustering. We also set Aq(u) := U(u), U n ■= U(u n ) and A q , n '■= A q (u n ), for all n £ N and 
q € No- Let 


_ h (Aq t n) 

»{U n ) ' 


( 7 ) 


Let B G B be an event. For some s > 0 and l > 0, we define: 


|_sJ+max{|/J —1, 0} 

W4B) = Pi r\B% ( 8 ) 

*=kl 


The notation f~ l is used for the preimage by f l . We will write W^ ( {B) := {yF s ^(B)) c . 
Whenever is clear or unimportant which event B £ B applies, we will drop the B and write 
just W s x or Wfp. Observe that 

Wo,n(U (u)) = {M n < u} and / _1 (^o,n(5)) = {^s > n}. 


(9) 
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Condition (/fq(un))- We say that JXfun) holds for the sequence Xq. X \,... if for every £,t,n £ 

N, 

n W tj £ (A q! n)) - H {A q , n ) H {%,£ (ig,n))| < 7ObM), (10) 

where 7 (q,n,t) is decreasing in t for each n and, there exists a sequence (tn)neN such that 
t n = o(n) and nj(q, n, t n ) —> 0 when n —> 00 . 


For some fixed q £ No, consider the sequence (t n ) n& fq, given by condition A(un.) and let 
(^n)neN be another sequence of integers such that 

k n —> 00 and k n t n = o(n). ( 11 ) 


Condition (ffi q (u n )). We say that fX q {u n ) holds for the sequence Xq, X\. Xq ,... if there exists 
a sequence (fcn)neN satisfying ( 11 ) and such that 


lim 

n—> 00 


\n/k n \ — l 

n ^2 h (A,n n f~ 3 {Aq,n)) = 0 - 

j=q+l 


( 12 ) 


From [FFT151 Corollary 2.4] follows that if the stochastic process satisfies both conditions 
JXq{u n ) and jX q (u n ) for some q £ N 0 , where (u n ) n£ n satisfies 0, then limn^oo fi{M n < u n ) = 
e~ 9r , whenever the limit 6 = hin^^oo 0 n exists. 

In this approach, it is rather important to observe the prominent role played by condition 
JX q {u n ). In particular, note that if condition /Xq( u n) holds for some particular q = qo £ No, 
then condition jX q (u n ) holds for all q > qo, which also implies that if the limit of 9 n in 0 
exists for such q = go it will also exist for all q > qo and takes always the same value. This 
suggests that in trying to find the existence of EVL, one should try the values q = qo until we 
find the smallest one that makes jXq{u n ) hold. With that purpose, as in |FFT151 Theorem 3.7], 
we consider the following. Let A £ B. We define a function that we refer to as first hitting 
time function to A, denoted by '■ X —> N U {+oo} where 

ta{x) = min {j £ N U {+oo} : f J (x) £ ^4} . (13) 

The restriction of 774 to A is called the first return time function to A. We define the first 
return time to A, which we denote by R(A), as the infimum of the return time function to A, 
i.e., 

R(A) = inf rA(x). (14) 

x£ A 

Assume that there exists q £ No such that 

q := min [j £ No : lim R(A^) = 00 } . (15) 

l n—Kx) J 

Then such q is the natural candidate to try to show the validity of /X q (u rl ). 


2.3. Rare Events Point Processes. A more sophisticated way of studying rare events con¬ 
sists in studying Rare Events Point Processes (REPP). These point processes keep record of 
the exceedances of the high thresholds u n by counting the number of such exceedances on a 
rescaled time interval. For every icMwe define 

^U n (A) := ^ ( ^-Xi>U n ■ 

ieAnNo 

In order to provide a proper framework of the problem we introduce next the necessary 
formalism to state the results regarding the convergence of point processes. We recommend 
the book of Kallenberg |Kal 86 | for further reading. 









CLUSTERING OF EXTREME EVENTS CREATED BY MULTIPLE CORRELATED MAXIMA 


7 


Consider the interval [0, oo) and its Borel cr-algebra £>[ 0oo ). Let x\,X 2 ,-- - £ [0,oo) and 
define 

OO 

v = ^2 s xi, 
i =1 

where 5 Xi is the Dirac measure at Xi, i.e., for every A G B \ 0)OO ), we have that 5 Xi (A) = 1 if 
Xi G A or 5 Xi (A) = 0, otherwise. The measure v is said to be a counting measure on [0, oo). Let 
A4 P ([0, oo)) be the space of counting measures on ([0, oo), £>[o l00 ))- We equip this space with 
the vague topology, i.e., v n —> v in _M p ([0,oo)) whenever u n (ijj) -A is(ip) for any continuous 
function if: : [0, oo) -A M with compact support. A point process N on [0, oo) is a random 
element of A4p([0, oo)), i.e., let (X,Bx, k) be a probability space, then any measurable map 
N : X —> A4 P ([0, oo)) is a point process on [0, oo). 

To give a concrete example of a point process, which in particular will appear as the limit 
of the REPP, we consider: 

Definition 2.2. Let T] . T 2 ,... be an i.i.d. sequence of r.v. with common exponential dis¬ 
tribution of mean 1/9. Let D\,D 2 , ... be another i.i.d. sequence of r.v., independent of the 
previous one, and with d.f. n. Given these sequences, for J G £>[ 0 , 00)1 set 

OO \ 

yy DiS Tl +...+Ti J, 

where 5t denotes the Dirac measure at t > 0. Let X denote the space of all possible realisations 
of Ti, T 2 ,... and D \, D 2 , ■ ■., equipped with a product cr-algebra and measure, then N : X —> 
A4p([0,oo)) is a point process which we call a compound Poisson process of intensity 9 and 
multiplicity d.f. ir. 

Remark 2.3. Throughout the paper, the multiplicity will always be integer valued which 
means that ir is completely defined by the values 7r f. = P(Di = k ), for every k G No- Note 
that, if 7Ti = 1 and 9=1. then N is the standard Poisson process and, for every t > 0, the 
random variable iV([0,f)) has a Poisson distribution of mean t. 

In order to define the REPP we need to rescale time. We do it using the factor v n : = 
l/P(Ao > u n ) given by Kac’s Theorem. However, before we give the definition, we need some 
more formalism. Let S denote the semi-ring of subsets of whose elements are intervals 
of the type [a, 6), for a, b G Mq\ Let TZ denote the ring generated by S. Recall that for 
every J £ TZ there are k G N and k intervals I\,... ,/& G S such that J = U^Lj Ij. In order 
to fix notation, let a,j,bj G Mq" be such that Ij = [dj,bj) G S. For I = [a, b) £ S and 
a £ M, we denote al := [aa,ab ) and I + a := [a + a, b + a). Similarly, for J £ TZ define 
aJ := al\ U • • • U al & and J + a := (Ii + a) U • • • U (/& + a). 

Definition 2.4. We define the rare event point process (REPP) by counting the number of 
exceedances during the (rescaled) time period v n J £ TZ, where J £ TZ. To be more precise, 
for every J G TZ, set 

N n {J) '■= • y/ Kl n { V n-I) = y ^ 1 Xj>u n - (16) 

i&v n JOHo 

We will see that the REPP considered here converge to a compound Poisson process. For 
completeness, we define here what we mean by convergence of point processes (see [Kal86] for 
more details). 
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Definition 2.5. Let (iV n ) ng pj : X -A _M p ([0, oo)) be a sequence of point processes defined on 
a probability space (X,Bx,d) and let N : Y — > _M P ([0, oo)) be another point process defined 
on a possibly distinct probability space (Y, By , v) ■ We say that N n converges in distribution 
to N if /j o iV " 1 converges weakly to u o iV _1 , i.e., for every continuous function (p defined on 
A4 P ([0, oo)), we have liirq^oc f pdfi o A ^" 1 = f tpdp o N _1 . We write N n =^=> N. 

Remark 2.6. It can be shown that ( N n ) n g p} converges in distribution to N if the sequence 
of vector r.v. (N n (Ji), ... , A r n ( ,/;,.)) converges in distribution to ..., N(Jk)), for every 

k G N and all J±,..., J\. £ 5 such that N(dJi) = 0 a.s., for i = 1,..., k. 


Note that 

{A„([0,n)) = 0 } = {M n <Un}, (17) 

hence the limit distribution of M n can be easily recovered from the convergence of the REPP. 


2.4. The convergence of REPP. In [ FFR j, the authors are preparing a synthesised version 
of the conditions in |FFT10| . in the absence of clustering, and in |FFT13j . in the presence 
of clustering, to a couple of general conditions that apply to general stationary stochastic 
processes, both in the presence and absence of clustering, which allow to prove the convergence 
of REPP. Moreover, these conditions, as the previous JJ^(u n ) and JX( u n), are particularly 
tailored to the application of dynamical systems and follow from a strong form of decay of 
correlations (against L 1 observables) to be defined below. 

Before we state the next condition we introduce some notation which will be useful in the 
rest of the paper. Let (l/l K - ) (u)) f . >Q be a sequence of nested open sets defined by: 

U (0 \u) := U[u) 

q 

4°)(u) := {X 0 >u,X i < u,...,X q <u} = U<®{u) Pi f| /“* (V°)(u)) c ) . 

i =1 

q 

U^ K \u) := U^-'Xu) - A^Xu) and MJ k) (u) := U^ K \u) P f] /"* ((f/ (re) (u)) c ) . 

2—1 


Let 


On := 




( 0 ) V 

n) 


l l ( LL 


and ir„ 


k = 


/x(4S) - /^(4?n +1) ) 


We notice that the set AqX was defined in Section 


2.2 


/^(4°n) 

and named A q , n . Since we need more 


Condition (^ 9 (tt n )*). We say that J\ q (u n )* holds for the sequence Xq,Xi,X 2 , ... if for any 
integers t, Ki,...,n and any J = U^ =2 /j G 1Z with inf{x : x G J} > t, 



where for each n we have that 7 (n,t) is nonincreasing in t and 717 ( 71 , t n ) —> 0 as —> 00 , for 

some sequence t n = o(n). 


In (FFRl Corollary 1 ] the authors show that under conditions J\ q (u n )* and /X q ( u n)-, where 
( u n ) n satisfies ([3|, the REPP (Al n ) ng ^ converges in distribution to a point process N whenever 
the limits 6 = lim n _^ 0O 0 n and 7 r(/c) = lim^^oo ir n (n) exist (for all k£N). 

If one assumes that the maximum of / occurs in a repelling periodic point £ and the 
probability ji is sufficiently regular then 
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(R2) the fact that ( is repelling means that we have backward contraction implying that 
there exists 0 < 9 < 1 so that P|* =0 /~ :ip (Xq > u ) is another ball of smaller radius 
around ( with 


h [ p| f~i p (X o > u) 
d=o 


(1 - 0)V(*o > «), 


for all u sufficiently large. 

Condition ( R2 ) guarantees that 9 n = q 'V converges to 9. Moreover, by |FFT131 Theo- 

H[Un) 

rem 1] the REPP ( N n ) n converges in distribution to a compound Poisson process N with 
intensity 9 and multiplicity distribution given by tt(k) = 9(1 — 9) K ~ l . 


3. Multiple correlated maxima 


In this paper the novelty of the approach resides inn the fact that instead of considering 
observables p : X —> M achieving a global maximum at a single point ( G X, we assume 
that the maximum is achieved at N points denoted by £i, £2 • • •, £/v> where in almost all cases 


(except in Sections 4.4 and 5.3), all these points are correlated by belonging to the orbit of a 
certain point £ € X, i.e., there exist mi < m 2 < ... < m n G Nq such that 


& = mo, 


(18) 


after possible relabelling of the maximal points. Let XI := {£ 1 , £2 ..., £at} be the set of 
maximal points of p. We will see that the binding between the maximal points expressed 


in (18) is responsible for a fake periodic behaviour which creates clustering of exceedances. 


Moreover, this mechanism of creating clustering turns out to be more flexible and allows to 
obtain easily different multiplicity distributions for the limit of the REPP (we recall that 
at repelling periodic points one obtains a Geometric multiplicity distribution) and different 
clustering patterns (we recall that at repelling periodic points the clusters consist of a strictly 
decreasing bulk of exceedances occurring at precise periodic intervals). 


3.1. Further assumptions on the observable. Besides having N global maxima £ 1 ,..., £jv, 
we assume that the observable p also satisfies 


p(x) = /ij(dist(x,£j)), \/x <E B £i (£i), i G {1, 


(19) 


where B Si (£i) n B e .(£j) = 0 and the function hi : [0,+ 00 ) —> M U { 00 } is such that 0 is a 
global maximum (/tj(0) may be + 00 ), hi is a strictly decreasing bijection hi : V —> W in a 
neighbourhood V of 0; and has one of the following three types of behaviour: 

(1) Type 1: there exists some strictly positive function g : W —> M such that for all y G M 


lim K x (s + yg(s )) 

s ~' > 9l(0) /lj (s) 


= e~y- 


( 20 ) 


(2) Type 2: hi(0) 


+00 and there exists j3 > 0 such that for all y > 0 


lim 

S—> OO 


hj \sy ) 

K H s ) 



(21) 
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(3) Type 3: hi(0) = D < +oo and there exists 7 > 0 such that for all y > 0 


lim h ‘) D ~ ^ = y\ 
h~ l {D-s) 


(22) 


Remark 3.1. Note that as long as the invariant measure has no atoms, then under the 
assumptions above on the observable we have that condition (Rl) is easily satisfied. 


3.2. Assumptions on the system and examples of application. We assume that the 
system admits a first return time induced map with decay of correlations against L 1 observ¬ 
ables. In order to clarify what is meant by the latter we define: 


Definition 3.2 (Decay of correlations). Let Ci,C 2 denote Banach spaces of real valued mea¬ 
surable functions defined on X. We denote the correlation of non-zero functions (j) G C\ and 
if G C 2 w.r.t. a measure P as 


Corp(0, if, n) := 


McAmc 2 


(j) (if o f n ) dP - / 0dP / if dP 


We say that we have decay of correlations, w.r.t. the measure P, for observables in C\ against 
observables in C 2 if, for every cf 6 C\ and every if E C 2 we have 


Corp(0, 1 /), n) —> 0, asn-> 00 . 

We say that we have decay of correlations against L 1 observables whenever this holds for 
C 2 = L 1 (P)and|^||c 2 = ||V>||i = /IV’|dP. 

If a system already has decay of correlations against L 1 observables, then by taking the 
whole set X as the base for the first return time induced map, which coincides with the 
original system, then the assumption we impose on the system is trivially satisfied. Examples 
of systems with such property include: 

• Uniformly expanding maps on the circle/interval (see |BG97| 1: 

• Markov maps (see ( BG97| ): 

• Piecewise expanding maps of the interval with countably many branches like Rychlik 

maps (see |Ryc83| ); _ 

• Higher dimensional piecewise expanding maps studied by Saussol in [SauOOj . 


Remark 3.3. In the first three examples above the Banach space Ci for the decay of correla¬ 
tions can be taken as the space of functions of bounded variation. In the fourth example the 
Banach space C\ is the space of functions with finite quasi-Holder norm studied in [SauOOj . 
We refer the readers to |BG971 ISauOOj or |AFV14] for precise definitions but mention that 
if J C M is an interval then 1 j is of bounded variation and its BV-norm is equal to 2, i.e., 
||lj|| w = 2 and if A denotes a ball or an annulus then 1a has a finite quasi-Holder norm. 


Although the examples above are all in some sense uniformly hyperbolic, we can consider 
non-uniformly hyperbolic systems, such as intermittent maps, which admit a ‘nice’ first return 
time induced map over some subset Y C X, called the base of the induced map. To be 
more precise, consider the usual original system as / : X —» X with an ergodic /-invariant 
probability measure fi, choose a subset Y C X and consider Fy : Y —> Y to be the first return 
map f ry to Y (note that F may be undefined at a zero Lebesgue measure set of points which 
do not return to Y, but most of these points are not important, so we will abuse notation here). 
Let //y(-) = be the conditional measure on Y. By Kac’s Theorem fj,y is Ty-invariant. 
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In |BSTV03] , the authors show that for typical points ( G Y the HTS to balls around 
( is the same both for the original and induced maps. This means, in particular, that the 
existence of EVL for the induced map implies the existence of an EVL for the original system, 
at typical points, and the limit is actually the same. In |FFT13j . a similar statement was 
obtained but this time for the convergence of the REPP at periodic points having the same 
compound Poisson limit. In |HWZ14| the authors show that regardless of the point ( G Y 
taken, the existence of a limit for the HTS (or EVL) for the induced first return time map 
implies the existence of a limit for the HTS (or EVL) of the original system, at f, and the 
limits coincide. This was generalised for the convergence of REPP in |FFTV15| . Namely, 
setting v% = 1 /py(X > u n ). for the induced process X} , 

£ R r>«„' 

j&vfJr iN 0 

Denote the speeded up return time by ta,y and the induced measure on Y by py. 

Theorem 3.4 ( |FFTV15| Theorem 3]). For r / > 0, setting J v := U s& jB v (s), we assume that 
N(J V ) is continuous in r], for all small r). 

N as n —> oo implies N n N as n ^ oo. 


This means that if we have a system that admits a first return time induced map with decay 
of correlations against L 1 , then as long as the points £i,... ,£j\r £ F, the convergence of the 
REPP (which implies the existence of EVL or HTS) for the original systems is determined 
by the convergence of the corresponding REPP for the induced map. Hence, we are reduced 
to proving the convergence of the REPP for systems with decay of correlations against L 1 
observables. This fact motivates the following: 


Assumption A Let f : X -A X be a systems with summable decay of correlations against 
L 1 observables, i.e., for all ip G C\ and if G L 1 , then Cor(ip,ip,n) < p n , with Yln>N Pn < oo. 
Moreover, we assume that there exists C > 0 such that for all n G N, we have 1 ,( 9 ), 1 .(«) £ C\ 

Fi-n <A-q,n 

and || 1 (,) Hcd || 1 .(«) ||ci < C. 

'S\q,n 


Among the examples of systems with these ‘nice’ induced maps we mention the Manneville- 
Pomeau (MP) map equipped with an absolutely continuous invariant probability measure. 
The form for such maps given in il.SVf)!). IBSTV03] is, for a G (0,1), 

, Jx (l + 2 a x a ) for xE [0,1/2) 

J MX) ~ [2x- 1 for x G [1/2,1] 

Members of this family of maps are often referred to as Liverani-Saussol-Vaienti maps since 
their actual equation was first introduced in |LSV99j . Let V be the renewal partition, that is 
the partition defined inductively by Z G V if Z = [1/2,1) or /(Z) G V. Now let Y G V and let 
Fy be the first return map to Y and py be the conditional measure on Y. It is well-known 
that (Y, Fy, py) is a Bernoulli map and hence, in particular, a Rychlik system (see |Ryc83| or 
[AFV141 Section 3.2.1] for the essential information about such systems). 
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3.3. Existence of distributional limits under Assumption A. As explained in |Frel31 
Section 5.1], conditions JX, q (u n ) and J\ q (u n )* are designed to follow easily from decay of correla¬ 
tions. In fact, if we choose (f = 1 _ 4 q n and = ly// Q e ^ q n j, in the case of Aq( u n), and <f> = 1^(«) 

and if = l n 4 :l )=n :j ■ i n the case of J\ q (u n )*, we have that we can take 7 (n,t) = Cpt- 

Hence, by taking (t n ) n such that linin-^oo npt n = 0, we get that conditions JX q {u n ) and Aq( u n)* 
are trivially satished. 

Now, we turn to condition JX q {u n ). Taking cb = fj = 1 ^ q n and since ||l^ 9i „||ci < C we 
easily get 

p(U n nf ■i(U n )') < (p(A q) n)) + || || ^-Aq t n |Il1(/z) Pj — + Cp(A.q, n )pj. 

(23) 


Let R n := R{A q) n)- Using estimate (23) and np{U n ) —> r as n —> 00 it follows that there 
exists some constant D > 0 such that 


[n/k n J L n/k n \ 

Tl ^2 n f~ j (-A q ,n)) = n ^2 MA,™ n f~ j (Aq >n )) 

j=q+ 1 j=Rn 

\n/k n \ / / » \\2 00 

^(Aq.n ) 2 + nCp(A q ,n) ^ Pj < ^ ^ - \~ U C p(Ag t n) ^ Pj 

j=Rn " U j=Rn 


if we check that linu^oo R n = + 00 . 


-> 0 , 


4. Maxima along a non-periodic orbit 


In this section we consider the case where ( is a non-periodic point. Our goal is to see what 
are the effects on the REPP when we compare it with the process obtained by an observable 
which attains its maximum at a unique point. As we will see, despite ( being non-periodic, 
the existence of multiple maxima along its orbit affects, for example, the extremal index and 
the multiplicity distribution. 

Note that by (-R1), we know that, for n sufficiently large, U n is the union of N disjoint balls 
around each that we denote by C/ n (£i) = (fit)- for i = 1,..., JV, so that 

N 

U n = \J UAfy). (24) 

i=l 


We begin by choosing a value of q for which condition fX{u n ) can be checked. From the 
computations in Section 3.3 we need to verify that R(A q , n ) —> 00 , as n —00 for such choice 
of q. 


Lemma 4.1. Assume that f satisfies Assumption A and that ( is a non-periodic point and 
that f is continuous in every point of the orbit of namely £, /(C), / 2 (C)--- Let An be 
as in (18). Let q = mjy — m\, then lim,^^ R(A qyn ) = 00 . 


Proof. Let q = m.jy — m\. By definition of A q: n we have that f J (-A q , n ) D A qt7l C f J {Aq,n) D 
U n = 0, for all j = 1,... ,q. Hence, we only need to check that the same holds true for all 
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j = q + 1,... ,q + L, where L is a given arbitrarily large integer. For that we will use the 
continuity of / over the orbit of £. Let 

Using the continuity over the orbit of £ and in particular of £jv, for every 1 < i < L, consider 
the balls such that f l {Bsfi^N)) C B e / 2 {f l {^N)) and consider the open set Dn = 

nf =1 B Si (£n)- Now, for each j = 1,..., N— 1 let Dj be a ball around £j so that f m N-m 3 (Dj) C 
D n . 

Finally, recalling (i?l) take n sufficiently large so that each U n (fi) C A n -®e/ 2 (£*) so 
that U n C U iL\Di n B e / 2 (£i)- Recall that A q<n C U n and then by construction for every 
j = 1,..., L we have that f q+q {A q: n) C U fLiB £ / 2 (f q+ ^(Ci)) and since by definition of e, we 
have B £ / 2 (f q+ ^(^i)) n U n = 0, then the result follows. □ 


4.1. The clustering effect and estimates for the extremal index. When the point £ is 
an absolute maximum and £ is non-periodic we know that the extremal index 8 is equal to 1. 
The situation is completely different when the maximum occurs at more than one point along 
the orbit of £. 


Proposition 4.2. Assume that f satisfies Assumption A. Let Xq,Xi,... be given by \ I), 
where the observable ip has N maximal points £i, ...,£jv related by (18) and satisfies (19), 
with £ E X a non-periodic point. In this case the extremal index 8 is given by 

Ell f (tf’«i) - Uli+1 


8 = lim 

n—>oo 


whenever the limit exists. 


z2i =i f 1 




Proof. In this case we have that q = mjsr — m\ and 




(o) _ 

q,n — 


where 


N 

^te) = c/ n 0) te)- U r {rnj ~ mt) (ui 0 ) m. 

j=i +1 


By the definition of 8 


lim 

n—)► oo 


h{Aq,l.) 

/i(c4 0) )’ 


the conclusion follows. 


□ 


In order to illustrate the richness of scenarios that the formula for the El of the previous 
proposition encompasses we consider below some particular examples. By condition (Rl), we 
know that, for n sufficiently large, U n = Un ^ is the union of N disjoint balls around each £, 
that we denote by Un' 1 (tfi). for i = 1,..., N, so that Un' 1 = U Depending on the 

different shapes of hi in ( |19[ ), the balls Un\Ci) may have different sizes and shapes. Moreover, 
if / is continuous along the orbit of £, then for 1 < i < j < N we have H 

Un (fij) 0 but, depending on the expansion rate and geometric distortion of / m >~ m ' 1 , it is 

not a priori clear if the one of the sets contains the other or not. Hence, in order to simplify, 
we assume that for 1 < i < j < N 

(i?3) C (uPHtj)) or (u^tij^ndu^tti) = 0, 

where the notation dB is used to denote the boundary of the set B. 
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Remark 4.3. Note that for all 1 < i < j < N, we have H Un\^i) / 

0. Hence, f-{ m j~ m i) n dUn\f.i) = 0 means that the connected component of 

f—(mj—mi) (Un\£j)J that intersects Un\^i) is completely contained in Un\^i). 

In the next result we treat the case where the condition (-R3) is valid. Before we assert the 
corollary we introduce some notation. 

Let I\ be the set of indices i E {1, ... IV — 1} such that 

N 

U n duW(ti) = 0. (25) 

e=i +1 


Under (R3), for every j 6 h, there exists ji E {i + 1,... N} such that 


N 


u r im - mi \u^m n c4°)(ei) = r (m ^- mi) (^ 0) (4)) n c4 0) (eo- 


e=i +1 


Corollary 4.4. If we are under the conditions of Proposition 4-2 and the condition (R3) is 
satisfied, then the extremal index 0 is given by 


= lim 

n—>oo 


Eiej, [(* («?’(&)) -M(/- (mj <-” ,) (C,f l (&))nf/< 0| (&))] + f .(ci 0, « K )) 


e«/*( cr<s 


r(°) / 


whenever the limit exists. 


Let us now assume we are in the particular case where /j is absolutely continuous with 
respect to the Lebesgue measure and its Radon-Nikodym density is sufficiently regular so that 
for all x E X we have 


lim 


p{B e (x)) 


dpi 


(x). 


(26) 


£->o Leb(LL(x)) dLeb' 

Note that if / is one dimensional smooth map and log(D/) is Holder continuous then, as seen 
in j FFTlhl Section 7.2], formula (26) holds. 


Corollary 4.5. If we are under the hypothesis of CorollaryUfjland p is absolutely continuous 


and (26) holds then the extremal index 9 is given by 



i =1 


Proof. We only need to notice that, by the mean value theorem and (26) 

1 




( r ^-™,) (f/ (o, tei)) n u L o, Kj) ) „ ______ 


□ 
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4.2. Estimating the multiplicity distribution. By Proposition 2.10 of |1FFR | we have that 
the multiplicity distribution of the limiting compound Poisson process for the REPP is given 
by the formula 


7T (k) 



(27) 


whenever the limit exists. 

Assume that (i?3) holds. 

Let /] be defined as in (25). For every i E we define o = *• Then, for 

i E {1,..., N — k} we inductively define If in the following way. 

Assuming that 4 is defined as well as jik-i, ■ ■ ■ > j/v-fc +2 fc-u then for every z = 1,..., N — 
k + 1, we define 4 & to the index such that 


AT 


U f (m ^" 1 mi) (^ 0) (W 1 )) n ^ 0) te) = / mi) (^ 0) fe,J)n^ 0) (6), (28) 


l=i +1 


when z E 4 and 4,fc = z, otherwise. 

Now, we define 4+1 as the set of indices z E {1,... IV — fc} such that 




U r (m ^ fc_m<) (c 4 0) (^, fc )) n= 0 . 

£=i+l 


Proposition 4.6. Assume that f satisfies Assumption A. Let Xq . X \, ...be given by Q, where 
the observable has N maximal points £i, ■.., £/v related by (|18|) and satisfies (19), with £ E X 


a non-periodic point. Assume that conditin (R3) holds. The multiplicity distribution of the 


point process is given by equation (27) (whenever the limit exists), where, for k < N, 


M4S) = a‘(/ {mjN ~ k ’ k mjv - fc) (c4 0) (^_ fc ,J)nc4 0) (e N -kj) 

t n c/fte)) - n n c/f(6)) 


+ E 

i&Ik+l 


and, for k > N, /i(Aq k l) = 0. 

Proof. We will prove by induction that, for i E 4+ 1 , 




4 fc+1) (&) = r imji ’ k+ ^ mi \u^^ jik+1 )) n 4 0) te), 

and, for z 0 4+i, /i (4?n(&)) = 0 since Un\&) = 










16 


D. AZEVEDO, A. C. M. FREITAS, J. M. FREITAS, AND F.B. RODRIGUES 


For i E I\ since j+o = £, we get 


N 


44&) = 4 0) (&)- U 


e=i +1 

IV 


f=i+l 

= Ofte) n „)) - r<*^M—>( 0 f(^ u )) > 

by definition of /,i. and, moreover, we have f i ' ib,; = / n f ,', " id, j 

N 

Note that, for i = N, = 0 and so 


£=?:+! 


4$(^) = ^ 0) (&v) = r {mj ^- mi \ui°\^ 0 )). 

For i£h,i±N, 4%) - USii/" (m ^ mi) (^n 0) te)) = 0, SO /X (4$(&)) = 0 and 

By induction, for i E 4+1 j we obtain 
N 

4Ste) = 4 fc) (&)- U r {me ~ mi) (u^m 

£=i +1 
N 

= ui k \ti) - U /-<"**-"*> (r im ^~ mi \ui 0 \^, k )) n t/fte)) , by induction 

£=i+i 

N 

= r {mji ^ mi \uL°\^ k )) n 4 0) te) - U /■ (m ^ _mi) (4 0) (^,J)> b y ^notion 

f=i+i 

= /- (m ^' m 44 0) (^,J) n 4 0) te) - by definition of j iifc+1 . 

and in particular we also have f7n fc+1) (4) = f~^ mji ‘ k+1 ~ mi \Un\^j ik+1 )) H Un\f,i)- 


N 


For i = N-k + 1, |J {U^ k) {^i)) = 0 and so 

l=i +1 


45(^-fc) = /' (m ^- fc+1 ’ fc_mAr - fe) (4 0) (^-, +1 ,J) n c4 0) (Civ-fc). 

For * 0 4+1, i ^ N - k + 1 p, (4^(6)) = 0 and 

4 fc+1) (&) = n 4 0) (6). 

For fc > iV + 1, /x (48(6)) = 0. 

This concludes the proof. □ 

Corollary 4.7. Under the assumptions of Proposition [7~h[ Assume that p zs absolutely con¬ 
tinuous with respect to the Lebesgue measure and (|26[) holds. Then the multiplicity distribution 
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of the point process is given by formula (27), where, for k < N, 

MMS) = Y ~ a n(m jik+1 ,mi)j +a n (m jN _ kk ,m N - k ) 


and 


*Tn(*L b'j — 


det Df a ~ b (f b (C))\ 


h 


Proof. Since p is absolutely continuous with respect to the Lebesgue measure and (26) holds, 
then 

1 




det Df mji ’ k ~ mi fa) 
Replacing this in the previous proposition we obtain the conclusion 


M (*4 0) (£*,*)) = a n (m ji k ,mi). 


□ 


Corollary 4.8. Assume that p is absolutely continuous with respect to the Lebesgue measure 
and (26) holds and 


f 


— (m i+ i—m;) 


(uf(£ m )) c uPiti 


Then the multiplicity distribution of the point process is given by formula (27) where, for 
k < N, 

m 43 )= E (“ n(m i+ k,mi) OLn (m^-j-fc+l; 777.^)^ T CXn(rflN:^N—k) 

*£-h +1 

and a n (a,b ) zs as m Corollary f.l 


Proof. Using the notation of Proposition 4.6, note that in this case ji k = * + k. 


□ 


4.3. An example. Let ( f,S 1 ,Leb,B ) be the dynamical system where / : S 1 —> S 1 is the map 
given by f(x) = 2x mod 1. Take ( = and an observable <p : S 1 —> M U { 00 } such that 

-log -(I, if x <E B € (C), 

I*-/(C) \Y ifx£B e (f(0) 

\x~f 3 ( C)|^, if x E B e (f(Q) 

„ 0, otherwie. 


¥>(z) = < 


where we choose e > 0, such that B e (f l (()) D B e (f J (()) = 0 for i,j £ {0,1,3} and i / j. In 
this case we have that uf = sup^ggi p(x) = 00 and it occurs at £, /(C) and / 3 (C). 

Given r > 0 and a sequence ( u n ) n satisfying ([3j we have that 

U n = (C - e-“», C + e"*") U (/(C) - u~\ /(C) + u~ 2 ) U (/ 3 (C) - u ~ 2 , / 3 (C) + a" 2 ). 

As before we let U n (C) = (C - e~“",C + e““”), t/ n (/(C)) = (/(C) - u~ 2 ,f(C,) + u~ 2 ) and 
U n (f 3 (0) = (/ 3 (C) — Ti/ 2 , / 3 (C) + w n 2 )- We notice that for n sufficiently large we have that 

(C - e~ Un , C + e~ Un ) g / _1 ((/(C) - w" 2 , /(C) + «n 2 )) , 

and 

r 2 ((/ 3 (C) - nc 2 , / 3 (C) + u~ 2 )) n t/ n (/(0) 5 (/(C) - ?v 2 , /(C) + u- 2 ). 

Noting that m\ = 0,m2 = 1,7773 = 3, then I\ = {2} and 72 = 3. So, by Corollary |4.7| we have 
that 
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8 = lim 

71—> OO 


I 1 1 

(4 0) (/ m2 ( o)) 

1 - n 1 

(r (m 3- m2 )(4 0 )(/ m 3(c)))) 

\ + p \ 

(4 0) (/ ms ( 0)) 

X/?;= i Z 2 1 

(4%)) 

1 


= lim 

n—Kx) 


' Dp(p(()) Un +U[ 

, -2 i —2 

e Un +U n + Un 



7 

8 ' 


To compute the multiplicity distribution we need to determine the sets Ik, for all k E N. It 
is easy to see that, in this case, I\ = {2} and Ik = 0, for all k > 2. We now need to determine 
the value of jj for all k E N and all i E /&. By definition j t Q = i, for all i E N. In addition, 
we can see that j' 2,1 = 3. Since mi = 0,m2 = l,m 3 = 3, by Corollary 4.7 we have 


M44 = («( ? «i 2 ,o,2) - at(m j2 ' 1>2 ) + a n (m j30 ,m 3 ))2u n 2 
= (a(m 2 , m 2 ) - a(m 3 , m 2 ) + a n (m 3 , m 3 ))2ri“ 2 

= 2 ^n 2 - Q) «n 2 + «n 2 

M44 = a (m j21 , m 2 ) = a(m 3 ,m 2 ) = 2u~ 2 , 

= a(m J12 ,mi) = a(mi,mi) = 2e~ Un , 

M44 = 0 for fc > 3. 

Therefore, 



7 T( 1 ) = lim M4%) = lim 


«n 2 - I A I “nHC-ld ^ 


—2 , 2 


-2 




( 0 ), 


—2 


—2 , „ —2 


- 1 A I Un + U n 


7t( 2) = lim M( 4 !n) ^44 = lim 

n ^°° /x(4in) n ^°° 

71(3) = lim ~ m( 4 S) = lim 

n—i►oo y(i(wA ^ ) n-^oo 

7r(fc) = 0 for A: > 4. 


2 J «n 2 -e u " 


«r ,. 2 - I A I +U 


—2 1 2 


1 

7’ 


1 “n 


«r ,. 2 - I A I UP+U. 


= 0 , 


—2 1 „ -2 


1 “n 


6 

7’ 


4.4. Non correlated maxima along non periodic orbits. Although the main purpose of 
the paper is to study the effect of multiple maxima when they are bound by belonging to the 
same orbit, we note that if we suppose that ip achieves a global maximum value at the points 
£ 1 ,... ,£ n but these maximal points have no intersecting orbits then one can also obtain a 
limiting EVL, which in the case of the maximal points being non-periodic, has no clustering 
and hence the El equals 1. 
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Let 0(x) denote the sequence of points that form the orbit of the point x, i.e., 0(x ) = 
{x,f(x),f 2 (x),...j . 

Proposition 4.9. Assume that f satisfies Assumption A. Let Xo,Xi,... be given by ([I]), where 
the observable (p has N maximal points £ 1 ,... , £jy Assume that all the maximal points are 
non-periodic, that f is continuous on U and that 

N 

n = 0 - ( 29 ) 

i= 1 


Then Xq,X\, ... satisfies conditions Jfo(u n ) andffi 0 (u n ) and consequently for (rt n )neN satisfying 
([3} then lim^oo /x(M n < u n ) = e~ r . 


Proof. In this case we take q = 0, meaning that Aq t n = U n . Conditions ,Ho( u n) and JXo( u n) 
follow from the computations in Section 3.3 as long as one check that linu^oo R{U n ) = 
+oo. This last statement follows from a continuity argument very similar to the proof of 
Lemma 14.11 □ 


Remark 4.10. This is in agreement with the result obtained in |HNT12| for multiple uncorre¬ 
lated maxima chosen as independent typical points of the invariant measure of non-uniformly 
hyperbolic systems admitting a Young tower with summable tails. Here, if the system is 
continuous and satisfies Assumption A (which is strictly contained in the class of systems 
considered in (HNT12| ). the statement can be reinforced since it applies to all non-periodic 
points chosen for maxima as long as they have non-intersecting orbits. 


5. Multiple maxima lying on periodic orbit 


In this section we consider the case where ( is a repelling periodic point of prime period p, 
meaning that there exists p 6 N such that f p (() = (, being that p is the smallest integer for 
which that happens, and the derivative Df p {C,) is defined and all its eigenvalues are strictly 
larger than 1. Note that in this case, if for example p is absolutely continuous with respect 


to Lebesgue measure and (26) holds then condition (R 2 ) is satisfied. From |FFT12| we know 
that if the global maximum of p is attained only at £ then we expect the existence of an El 
equals to 

1 


0 = 1 - 


det DfP(Q \ ’ 


in the case p is absolutely continuous with respect to Lebesgue and (26) holds. Moreover, by 
|FFT1 3] . we know that the REPP converges to a compound Poisson process with a geometric 
multiplicity distribution given by 


tt ( k ) = 6{i-ey~ l , 

for all k > 1. 

In this section, our goal is to study the effect of having multiple maxima along the orbit 
of the periodic point £. We will see the El is affected (it decreases) and asymptotics of the 
REPP is also affected. 

We assume without loss of generality that tp achieves a maximum at the points £ 1 ,... An, 
where each is given by equation (18), where 0 = mi < m 2 < ... < mjv < p — 1, so that, in 
particular £1 = (. 
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We begin by choosing a value of q for which condition JX(u n ) can be checked. One easily 


anticipates that a suitable choice is taking q = p. From the computations in Section 3.3 
need to verify that R(A qt n) -A oo, as n —> oo for such choice of q. 


we 


Lemma 5.1. Assume that £ is a repelling periodic point of prime period p. Let £i,... ,£jv be 
as in (18), where 0 = m± < m 2 < ... < m n < p — 1. Let q = p, then liirp^oc R(A q , n ) = 00 . 


Proof. Let g := f p , and so every is a fixed point of g. As ( is a repelling periodic point, 
we have that there is a neighbourhood of fj that we denote by V(fi) where g |\/(^) : P(£i) —> 
g{V{ii)) C X is a diffeomorphim. By the Hartman-Grobman theorem there is a neighbour¬ 
hood W{ 0) of 0 E T^X and a homeomorphism h : W(0) —> P(£) such that h o Dg((f) = g oh. 

As W(ff) is an open set, there is r > 0 so that the ball of radius r, denoted by B r { 0) 
is contained in W{ 0). Let x E B r (£i) and let e = dist(x, £j). By hypothesis, the point 
is repelling periodic and it implies that there exists A > 1 such that the diameter of 
(Dg(()) K (B e (£ i i)) is at least A K e. So, in order for X K > r , we must have that n > log iogA S£ - 
Hence, Dg(())«(B e (&) C B r (&), for « < 

Now, for each j E {1,, iV}, let V(£j) be an open neighbourhood of where f\v (£■) : 
V(fj ) —>■ V(f mj+1 (C)) is a diffeomorphism and the Hartman-Grobman theorem holds. Let n E 
N such that U n (fj) C V (f-j ) . Given x E Aq :ni by the definition of A q>n we know that g K (x) fi 
Aq t n while g K (x) E B r (fi). Let e n be the diameter of A q . n . Then while k < Iog £n we have 

that g K (x ) ^ A q j,,. Since log£ ” -+ 00 , as n —> 00 , then it follows that lin^^oo R(A qjn ) = 
00 . □ 


5.1. The extremal index. In what follows we need some notation, which we introduce now: 
We define mi + jv := mi + p, so that f mi+N (() = f mi { C). 


Proposition 5.2. Let Xq. X\,... be given by ([!]), where the observable is given by (19), with 
a repelling periodic point. In this case the extremal index 6 is given by 


N 


N+i 


= lim 

n—> 00 


whenever the limit exists. 


j=i+l 


2=1 


N 


E" 

2=1 


Proof. Note that A^} n = U^li where 

N+i 

4S&) = - U 

j=i +1 


Since 6 = lim 

n—^ 00 




the conclusion follows. 


□ 


Just like in the non-peridic case, we now consider some particular examples that illustrate 
the richness of scenarios that the formula for the El of the previous proposition encompasses. 
Hence, as before, we assume that condition (R3) holds. 

First, we introduce some notation. 
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Let I\ be the set of indices i G {1,... N} such that 

N+i 

U n du^^i) = 0. 

t=i +1 

Under (-R3), for every '(£/]. there exists ji G {i + 1,... N + i\ such that 
N+i 

U n c4°>(£i) = /-("*«-^ ) ( c4°)(^.)) n 

2=i+l 


(30) 


Corollary 5.3. If we are under the conditions of Proposition 5.2 and the condition ( R3 ) is 
satisfied, then the extremal index 9 is given by 


6 = lim 

n—>oo 


Sie/i 


{U [ n\fi)) 

~ P 

(C 






1 p 




whenever the limit exists. 

Let us now assume we are in the particular case where p is absolutely continuous with 
respect to the Lebesgue measure and (26) holds. 


Corollary 5.4. If we are under the hypotheses of Corollary 5.3 and p is absolutely continuous 

1 


and (26) holds then the extremal index 6 is given by 

P 


9 = lim 

n—toc 


E 

ieh . 


I det £>/"*«-"*(&) I 


P 


n \S,Ji 


N 

Em 

2=1 

Proof. The proof is analogous to the proof of the non-periodic case presented in Corollary |4.5[ 

□ 

5.2. The multiplicity distribution. As mentioned before, by Proposition 2.10 of |FFR| we 
have that the multiplicity distribution of the limiting compound Poisson process for the REPP 
is given by the formula (27), whenever the limit exists. 

Recall that we are assuming that (-R3) holds. 

Let I\ be defined as in (30). For every i G {1,...,1V}, we define = i. Then, for 
i G {1 ,... ,N} we inductively define Ik in the following way. 

Assuming that Ik is defined as well as ji k—i, ■ ■ ■, jN,k-i • then for every i = 1,..., N, we 
define ji+ to the index such that 
N+i 

U n ^ 0) fe) = <) ( t4°)(e iiifc)) n (31) 

e = i +1 

when i G Ik and ji k = i, otherwise. 

Now, we define Ik+i as the set of indices i G {1,... N — k} such that 
N+i 

U J) n duPiti) = 0. 

2 = 1+1 
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Proposition 5.5. The multiplicity distribution of the point process is given by eqiLation (??) 
(whenever the limit exists), where, for k < N, 


M (AS)= £ 

(-1 







and, for k > N, p,(Aq k }i) = 0. 

Proof. We will prove by induction that, for i £ Ik+ 1 , 

4$(&) = u^Xti) n r^-^Xut°x^. k )) - r {mj ^- mi Xui°x^. k+1 )) 

and 

4 fc+1) (£0 = r {mji ^~ mi Xui°X^ k+1 )) n u(°Xa 

and, for i 0 I k+1 , p = 0 since UnXfi) = Un^XCi). 

For i £ 1 1 since jig = £, we get 

N+i 

A< 0 l(t,) = U< 0 >(( i )- U 

l=i+1 
N+i 

= d, 0> ({.)- U r <m ».«- m - l (t 4 " ) ({* >0 )) 

l=i+ 1 


by definition of and, moreover, we have (£j) = / t )) n Un’^i). 


n\tii) = /- (m *,i-^>(c4 0) (^J) n t4 0) ( 

For i<f h,i A N, 4%) - UStr /- (w *«-”*i ) (c4 0) (Ci)) = 0, so // (4$(&)) = 0 and 

m( 4 1} (6)) =h(4 0) (e,))- 


By induction, for i £ dfc+i, we obtain 
Af+i 

4S&) = 4 fe) (&)- U 

i=i+1 
Af+i 


4 fc) (6.) - U r (m/ “ mi) (/- (m ^- m|) (c4 0) (^ ifc )) n 4 0} (^)) , by induction 


€=i+l 


- /" Ki ' fe ' m °(4 0) (4, fc )) n 4 0) te) -h( U r (m «.*" roi) (c4 0) (^, fc )) ) , by induction 


Af+i 


V.&=*+1 


= r {mji ’ k ~ mi) (Ui°X^. J) n 4 0) te) - /- (m W-*)(^ 0) ( U+i )) > by definition of 
and in particular we also have 4 fc+1) (&) = f~^ mji ’ k+1 ~ nH XUnX^j i jfc+1 )) H UnX^i)- 
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By induction, for i G I k +i, we obtain 


N+i 


e=i +1 


= c4 fc) (&) - U (/ _(m ^ fc - mi) (^ 0) (^,J) n t4 0) (&)) , by induction 

fc=i+l 

/ IV+* \ 

= r {m ^ mi \u^\^ k )) n U^\£i) - pt MJ , by induction 

\ e = i +i / 

= n - /-^A + i- m ^(C/( 0) (^. fc+1 )), by definition of j ijk+1 . 


and in particular we also have Un +1 \&) = f ( ' mji ’ k + 1 m '\Un\£j i:k+1 )) H Un’iZi )• 

For i 0 I fc+ i, i^N-k + ln = 0 and 

u£ +l \ti) = /- (mj ^- mi) (^ 0) (4,J) n e4 0) &). 

This concludes the proof. □ 

Corollary 5.6. Assume that /x is absolutely continuous with respect to the Lebesgue measure 
and (26) holds. Then the multiplicity distribution of the point process is given by formula (27), 
where, for k < N, 

M-4S) = S -«nK ii+1) m,)) +a n (m jN _ kk ,m N - k ) 

*£bc+1 


r(°) / 


and 




det D f a ~ b (f b (())\ 


l l 


Proof. The proof is analogous to the proof of the non-periodic case presented in Corollary |4.7| 

□ 


Corollary 5.7. Assume that /x is absolutely continuous with respect to the Lebesgue measure 
and (26) holds and 


f 


— (m i+ i—m;) 


(uL 0) (ti+ 1)) C uW(ti 


Then the multiplicity distribution of the point process is given by formula (27) where, for 
k < N, 

= ^2 (a n (m i+k , mi) - a n (m i+k+1 ,mi )) + a n (m N , m N _ k ) 

i&Ik+l 


and a n (a,b ) is as in Corollary 5.6- 


Proof. Using the notation of Proposition 5.5, note that in this case Ji k = i + k. 


□ 
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5.2.1. Example. Let ( f,S 1 ,Leb,B ) be the dynamical system where / : S 1 —> S 1 is the map 
given by f{x) = 2x mod 1. Take the periodic point of period 5, C = and an observable 
p : S 1 —> R U {oo} such that 


p(x) 


' -log|x-C|, 

l»- /(C) 1^, 

\x- / 3 (C)I^, 

0, 


if X e B e ( C), 

if x€B e (f(Q) 
ifxeB e (f(0) 

otherwie, 


where we choose e > 0, such that B f (f t (Q) PI B e (f J (()) = 0 for i,j G {0,1, 3} and i / j. In 
this case we have that uf = sup x . e gi p{x) = oo and it occurs in C, /(C) and / 3 (C). 

Given r > 0 and a sequence (iv)neN satisfying [3] we have that 

Un = ((- e~ Un , C + e _Un ) U (/(C) - I/; 2 , /(C) + «n 2 ) U (/ 3 (C) - / 3 (C) + U~ 2 ). 


We notice that for t sufficiently large we have that 

(C - e~ Un , C + e-“") 5 /- 1 ((/(C) - V} 2 , /(C) + n- 2 )) , 

/- 4 (C - e"“", C + e~ Un ) £ /" 2 ((/ 3 (C) - u~ 2 , / 3 (C) + IV 2 )) £ (/(C) - IV 2 , /(C) + IV 2 ) , 

and 

r 2 (c - e_Un , c+e-“”) s /' 3 ((/(O - <v 2 , /(o+iv 2 )) s (/ 3 (o - iv 2 , / 3 (0+ u ~ 2 ). 

The extremal index: Noting that mi = 0, m 2 = l,m 3 = 3, then Iq = {2,3} and 72 = 


3, J 3 = 5. So, by Corollary |5.6|we have that 


6 = lim 

n—> oo 


/ 

V 

+ 


h (t4 0) &)) - 


Df m 3~ m 2 (f m 3 (C)) 


/I 


m 3 


E£im(V>({,)) 


m 5 


E£i**(V”te)) 


= lim 

n—>■ 00 


u -2 _ 1 _ 

n Dp(p( 0) 


v 2 + \ 2 - 


DP(P(0) 


u. 


-2 


- 2,-2 


2-j-g 13 

2 16' 


e + U n +U~n 

The multiplicity distribution: To compute the multiplicity we will use Corollary |5.7[ 


First, we need to determine the sets p , for all k G N. It is not difficult to see that, in this 
case, p = {2, 3}, for all k G N. We now need to determine the value of jip . for all k G N and 
all i G Ik - By definition j/0 = / for all i G N. We can verify by induction that 

j'2.2t. = 3£ — 2, j2,2Vi = 3£ — 3, 

J3,2f = 3£ — 3, J3,2£—1 =3^ — 5. 

















CLUSTERING OF EXTREME EVENTS CREATED BY MULTIPLE CORRELATED MAXIMA 


25 


/m 

In order to compute *4.^, note that in this case q = 5. So, 

®Gifc+i 

= o:(m htk ,m 2 ) - a{m h k+1 ,m 2 ) + a(m j3k ,m 3 ) - a(m hh ,m 3 ). 
Considering first the case A; = 2Z, with f £ R, we have 

M-4 2 ?) = ot{m j22V m 2 ) - a(m J - 2i!M _ 1 ,m 2 ) + a(m j3 2e ,m 3 ) - a(m h 2£ _ 15 m 3 ) 

= a(m 3 £_ 2 , m 2 ) - a(m 3 £_ 3 , m 2 ) + a(m 3 ^_ 3 , m 3 ) - a(m 3 £_ 5 , m 3 ) 



For the case k = 2£ + 1 with / 6 No, we have 

MA 5 2 n +1) ) = a ( m h, 2 t- 1 ^ 2 ) - a(m h 2e _ 2 ,m 2 ) + a(m hi!U _ t , m 3 ) - a{m j3 2e _ 2 ,m 3 ) 
= a(m 3£ - 3 , m 2 ) - a(m 3 ^_ 5 , ra 2 ) + a(m 3 £_ 5 , m 3 ) - a(m 3£ - e ,m 3 ) 


5£-2 


^n 2 - I - 


5£-5 


+ 


5£-3 


^n 2 - I - 


5£-5 


= 


5 /I 


16 V 2 


5£ 


-2 


Replacing these in the expression in Corollary 5.6 we obtain 

U r / 1 \ 5£ 


M4“ ) )-M4;r i ')_TV2 


7r(2/+ 1) = 


13 /IV 1 _ 2 5/1 _ 

(2^+14 "7T ( 4 ) U n Z - — [- U r 


16 V 2 




13 


-u. 


-2 


21 /I 

26 \ 2 


u 


, for t G No, 


M4V 1 ’) - m( 4“>) _ i6 VI 


7t(2/) = 


5 ^lV (£_1) „.-2 13 ( l 


U n 8 l 2 ' Un 


5 1 


-2 


67 /I 


m45) 


13 


-u. 


-2 


13 V 2 


5 1 


, for le N. 


5.3. Non correlated multiple maxima including periodic points. We discuss briefly 


the possibility of having multiple maximal points with non-intersecting orbits as in Section 4.4 


but now allowing the possibility of having maximal points that are periodic. By [FFT131 The¬ 
orem 2] and |AFV14l Theorem A], we have that for every continuous system under Assumption 
A, for every £ € X taken as the only maximum of 93 , then there exists an EVL with an El 
that is 1 if the point is not periodic or is less than 1 (whose value depends on the expansion 
rate at £). Hence, we are going to consider multiple maximal points £ 1 , ■ ■ • ,£/v and to each 
such point we associate the respective El guaranteed by |FFT13l Theorem 2] and |AFV141 
Theorem A] (see also |Frel31 Theorem 8 ]), which we denote by 9 t . being that if £,; is periodic 
then 9j < 1 . 


Proposition 5.8. Assume that f is continuous and satisfies Assumption A. Let Xq,X\,... 
be given by Q, where the observable ip has N maximal points £1 , ...,£jv an d let 9% be the 
corresponding EL Let p % be the period of £* (we set Pi = 0 when £j is not periodic). Assume 
that (29) holds and that the following limits exist for every i = 1,..., N 


lim 

n—>00 


h{U n ) 
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where U n (fi) is a connected component of U n as in (24). Then Xq. X\,... satisfies conditions 
ffq( u n) and JX q {u n ), for q = maxj = i an d consequently for ('U n )neN satisfying ([3j then 
linin^oc l-i(M n < u n ) = e~ dr , where 6 = Y,iLi a ^i- 


Proof. Conditions Jf q {u n ) and ffi (u n ) follow from the computations in Section 3.3 as long as 


one check that linin^oo R(A qq i) = +oo. This last statement follows from the definition of A, 


and a continuity argument very similar to the proof of Lemma 4.1 


q,n 

□ 


6. Clustering patterns 

As we have seen in Sections [4] and [5j multiple correlated maxima can be used as a mechanism 
to create clustering, alter the value of El and produce different multiplicity distributions for the 
limit of REPP. When compared to the usual method of introducing clustering in a dynamical 
setting, which was essentially based on considering a global maximum of p achieved at a 
unique repelling periodic point, this method of forging clustering by means of considering 
multiple correlated maxima is much richer. 

In fact, we have already seen that, on the contrary to single maximum at a repelling periodic 
point for which one always obtains a geometric distribution for the multiplicity distribution 
of the limiting compound Poisson process for the REPP, with multiple correlated maxima we 
can obtain different multiplicity distributions for the limit of the REPP. 

But the richness of this mechanism can also be appreciated by simply looking at the observed 
data and the corresponding clustering patterns. Observe that the systems we are studying 
are either uniformly expanding or non-uniformly expanding, which means that the periodic 
points should be repelling. When the maximum is achieved at a single repelling point the 
clustering pattern is very simple: it consists of a very large observation corresponding to the 
first exceedance, which is followed by a strictly decreasing cluster of exceedances observed after 
a precise number of observations corresponding to the period (see Figure[l]). This is because in 
order to have a cluster of exceedances the orbit needs to enter in a very small neighbourhood 
of the periodic point ( of period p. Let us assume we have a very high exceedance, which of 
course means that the orbit went very close to f. Once inside this very small neighbourhood 
of (, the periodicity of ( forces the appearance of another exceedance after p observations but 
since ( is repelling this return to a small neighbourhood of ( is not that deep, i.e., the orbit 
is pulled away from f which explains the fading of the exceedances. 


FIGURE 1. Cluster of observations for an observable with global maximum 
achieved at the fixed point 1 /2 of the system / given by (32) 
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In order to illustrate how easily one can create a different clustering pattern, where in 
particular one can observe that the exceedances may grow inside a cluster, we consider the 
following example. 

We consider the system uniformly expanding map 

/ : S 1 — > S 1 (32) 

x i-a 3x mod 1 

and the observable ip : S 1 —> M, where <p(x) = — log|x — ||lj 0 ij + \x — || _ 3l^i 1 j. This 
observable the maximum value is uf = sup, cg gi p(x) = oo, which is achieved at the points 
£ = ^ and /(£) = |, being that £ = / 2 (£), ie., , the maxima correspond to the points of a 
periodic orbit of period 2. 

Note that for u sufficiently large, {Xo > u} = (1/4 — e~ u , l/4+e _,u )U (3/4 — w^ 3 , 3/4+u -3 ). 
So an exceedance of u occurs when the orbit enters a ball of radius e~ u around 1/4 or a ball 
of radius u~ 3 around 3/4. Note that if we have an exceedance because the orbit hits the point 
x = 1/4 + e -u /3 and p(x) = u + log 3 > u, then in the next iterate the orbit hits the point 
f(x) = 3/4 + e~ u and ip(f(x)) = e u / 3 u + log 3 > u, for u sufficiently large so that we have 
an exceedance that is followed by another exceedance that is even higher than the first. 

This way, we create different clustering patterns which can be fully appreciated by the 
simulation study we performed. We chose a random point x £ [0,1] (using the uniform 
distribution); we iterated the point 2000 times by /; we considered u to be such that e -u + ^ 3 - = 
(meaning that the expected number of exceedances of u is 40). In Figure [ 2 ] we can see the 
actual time series and in Figure [3] we can see the observations forming a cluster that presents 
clearly a different pattern, with exceedances being followed by even higher exceedances, which 
was forbidden in the case of a single maximum as portrayed in Figure [l] 





Figure 2. Time series for 2000 runs of the orbit of the point x 
0.7756592465669858 with the level u = 4.619613119957849 depicted in red 


7. Competition between domains of attraction 

Up until now, we have addressed the issues regarding the effect of having multiple maxima 
on the existence of EVL, the appearance of clustering, the consequences on the value of the El 
and the impact on the limit of REPP. This has been done for normalising sequences (u n ) ne pj 
such that condition ([ 3 ]) holds. Typically, these normalising sequences are taken as a one 
parameter family of linear sequences such that u n = y/a n + b n , where (a n ) nS N, (b n ) n gpj, with 
a n > 0, and one usually looks for non-degenerate limit distributions for P(a n (M n — b n ) < y ). 

The limit is then e~ T ^ y \ where r{y) is a function of y, which is connected to the type of 


observable we have. For well behaved measures, if the observable satisfies (20) then one gets 


that r(y) = e v meaning that the the stochastic process falls in the domain of attraction of 
a Gumbel law (or type 1); if the observable satisfies (21) then one gets that r(y) = y~ a (for 


y > 0 and a > 0) meaning that the the stochastic process falls in the domain of attraction 
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FIGURE 3. Cluster occurring between the observations 1446 and 1450 of the 
time series in Figure [2] 


of a Frechet law (or type 2); observable satisfies (22) then one gets that r(y) = {—y) a (for 
y > 0 and a > 0) meaning that the the stochastic process falls in the domain of attraction of 
a Weibull law (or type 3). 

In the independent and identically distributed (i.i.d.) setting the domain of attraction is 
determined by the tail of the distribution F. where F(v) = ^(Xq < u). As it was observed 
in jFFTlOl. Remark 1], in the dynamical setting (with a single maximum), the shape of the 
observable determines the type of tail of the distribution F and, once the existence is estab¬ 
lished by the assumptions on the system, then the domain of attraction is determined by the 
behaviour of the observable at the maximum. 

In here, since we have multiple maxima, then we may have different types of behaviour of 
the observable at the different maximal points. This creates a sort of competition between 
different domains of attraction, which in fact is not that difficult to settle. First, we start 
by noting that if the observable is of type 2 (satisfies (|2T|) ) then supt/? = h(0 ) = oo; if the 
observable is of type 3 (satisfies (|22[) ) then sup</? = h{ 0) = D < oo and if the observable is of 
type 1 (satisfies (20)) then the maximum value can be either finite or infinite. Hence, since 
all maxima achieve the same maximum value, we can only have competitions between Frechet 
and Gumbel and between Weibull and Gumbel. 

Note that the appearance of clustering and consequently of an El less than 1 does not 
change the type of limit law that applies when we consider one parameter linear sequences. 
This is the statement of [LLR831 Corollary 3.7.3], which also adds that, in fact, we can even 
perform a linear change of the normalising constants a n and b n so that the exact same limit 
applies as in the corresponding independent case. Hence in what follows we restrict the study 
to the behaviour of y,{U{u)) = 1 — F{u ) instead of /i(A^(u)) as u goes to supyx 

In what follows we establish the reigning domain of attraction in each such situation. For 
that purpose, we introduce the following notation F(x) = 1 — F(x). 


7.1. Gumbel versus Frechet. 
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Proposition 7.1. Let ip : X — > M U {00} be an observable which has the following form 

<p(x) = <pi(x)1u(q(x) + ip 2 {x)l u{fi ^(x), 

where <p\ is an observable defined in the neighborhood U(Q) of(, ip 2 is an observable defined in 
the neighborhood U(f' L ((f)) of /*(£) and they satisfy max x <pi(x) = ipi(C) = 00 = max x ip 2 (x) = 
<^2 (/* (C)) • -hef Fi be the distribution function of ipi. If F\ and F 2 belong to the domain 
of attraction for maxima of the Gumbel and the Frechet distribution, respectively, then the 
distribution function of ip belongs to the domain of attraction of the Frechet distribution. 


Proof. We start by noting that if F\ belongs to the domain of attraction of the Gumbel 
distribution, with infinite right endpoint of the support, then we may write 

1 - Fi(u) = u p e~ uq Li(u ), (33) 


where p E M, q > 0 and L\ is a function of slow variation at infinity, i.e., 
( 0 , +00) and, for all t > 0 , 


lim 

U —^OO 


L\{tu) 

L\(y) 


= 1 


L\ : (0, +00) — y 


(cf. |Gom84j ). 

F 2 belongs to the domain of attraction of the Frechet distribution if and only if 


1 - F 2 (u) = u a L 2 {u), 


(34) 


where a > 0 and L 2 is a function of slow variation at infinity (cf. IdHTOl h 
If 

tp(x) = (Pi(x)1u(q(x) + <p 2 (x)l U ( f i(Q' > (x), 

and Xq = <p, then 

fi(ip > u) = p,(X 0 > u ) = u p e~ u<1 L\{u) + u~ a L 2 (u ) = u~ a ( L 2 (u ) + u a+p e~ u<1 L\(u)) . 


As lim^oo u a+p e~ uq L\(u) = 0 we get that the function L(u ) = L 2 (u) + u a+p e~ u<1 L\(u) is a 
function of slow variation at infinity. Hence, the distribution function of <p(x) belongs to the 
domain of attraction of the Frechet distribution. □ 


7.2. Gumbel versus Weibull. 

Proposition 7.2. Let ip : X —> M U {00} be an observable which has the following form 

ip{x) = ¥>i(®)ltf (c) (a;) + cp 2 (x)l u{fi{0) (x), 

where ip\ is an observable defined in the neighborhood U(() of Q, ip 2 is an observable defined in 
the neighborhood U (/*(C)) °f /*(C) an d they satisfy max x p>i{x) = ip\{Cf) = D = max x ip 2 (x) = 
<P2(P(0)- Let Fi be the distribution function of ipi. If F\ and F 2 belong to the domain 
of attraction for maxima of the Gumbel and the Weibull distribution, respectively, then the 
distribution function of ip, F, belongs to the domain of attraction of the Weibull distribution. 

Proof. If Xq = ip, then 

1 - F(u) = n(ip > u) = p(X 0 > u) = (1 - Fi(u)) + (1 - F 2 {u)). 
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So, 


1 -F(D-sy)_ 1 - Fi(D - sy) + 1 - F 2 (D - sy) 

.s™ 1 -F(D-s) ~sTo 1 - F^D - s) + 1 - F 2 (D - s) 

(1 - F 2 (D - sy)) ( ilgjglgj + l) 


= lim 

s —>o 


(i - f 2 ( d -*)) ( !:g|g:;j + i 

As F 2 belongs to the domain of attraction of the Weibull distribution, then 

1 - F 2 (D - sy)) 

An Fjn si 

for some a > 0 , and 

1 - F 2 (D - s)) = s a L(s), 
where L is a function of slow variation at zero, i.e., for all t > 0, 

lim hid = 1. 

S —>0 LySj 

We note now that if F\ belongs to the domain of attraction of the Gumbel distribution, 
with finite right endpoint of the support D, then, for all (3 > 0, 

, l-FUD-s) 
lim- 3 - = 0 . 

s-i >0 sP 


(35) 

(36) 

(37) 


Thus, by (37), we have that 


, 1-FUD-s) 

hm- — — -- = 0 . 

s->o 1 — F 2 (D — s ) 


(38) 


Using now (36) and (37) in (35), we obtain that 


1 - F(D - sy) 

hm-——-— = y , 

s— i-o 1 — F(D — s ) 

and so F belongs to the domain of attraction of the Weibull distribution. 


□ 
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